Some Systems of Multivariable Orthogonal 
Askey- Wilson Polynomials 

George Gasper* and Mizan Rahmanf 
Abstract 

In 1991 Tratnik derived two systems of multivariable orthogonal Wilson polynomials 
and considered their limit cases, q- Analogues of these systems are derived, yielding systems 
of multivariable orthogonal Askey- Wilson polynomials and their special and limit cases. 

1. Introduction. In [10] Tratnik extended the Wilson [12] polynomials 

(1.1) w n (x; a, 6, c, d) = (a + b) n (a + c) n (a + d) n 

-n, n + a + b + c + d — 1, a + ix,a — ix 



X4F3 



;i 

a + o,a + c,a + d 



to the multivariable Wilson polynomials (in a different notation) 
(1.2) 

W n (x) = W n (x; a, 6, c, d, a 2 , 03, . . . , a s ) 

"s-1 

Y\w nk (x k ;a + a 2 ,k + N k -i,b + a 2 ,k + Nk-i,a,k+i +ix k +i,a k +i - ix k+1 ) 

.k=l 

x w na (x s ; a + a 2 , s + N 8 -i, b + a 2 , s + N 8 -i, c, d), 
where, as elsewhere, 

k 

(1.3) x= (xi,...,x s ), n = (m, . . . ,n s ), a jjk =/ J a >i, OL k =a x , k , 

k 

N o,k = Nk = Nl ' k > a k+i,k = N k+ljk = 0, 1 < j < k < s. 

These polynomials are of total degree iV s in the variables yx, . . . , y s with y k = x k , k = 
1, 2, . . . , s, and they form a complete set for polynomials in these variables. 

In Askey and Wilson [1], [2] the notations W n (x 2 ; a, 6, c, d) and p n (—x 2 ) are used for 
the polynomials in (1.1) and their orthogonality relation is given. Tratnik [10, (2.5)] proved 
that the VF n (x) polynomials satisfy the orthogonality relation 

/oo ^00 
• • • / W n (x) W / " m (x)p(x) dx X ■ ■ ■ dx s = A n 5 njln 
-00 J — 00 
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for Re(a, 6, c, d, a 2 , . 
(1.5) p(x) 



, a s ) > with 

- T(a + ixi)r(a — + ixi)T(b — ix\) 



x 



s-l 

n 

.fe=i 



r(a fc+ i + ixk+i + ixfc)r(a fc+ i - zxfc+i - ix k ) 
T(2ix k ) 



x 



X 



r(afc+i + ix fc+ i - zx fc )r(a fc+ i - ix k +i + ixfc) 
r(-2za: fc ) 

T(c + ix s )r(c — zx s )r(<i + ix s )T(d — ix s ) 
T{2ix s )T(-2ix s ) ' 



(1.6) A n = (4n)< 



s 

Y[ n k l (N k + N k _ x + 2a k+1 - l) nfc 
.k=i 



r(iV fc + N k ^ + 2a k )T(n k + 2a fc+1 ) 
T(2iV fc + 2a k+l ) 

x T(a + c + a 2)S + iV s )r(a + d + a 2 , s + N s )T(b + c + a 2 , s + N s ) 
x T(b + d + a 2 , s + N a ), 

and 2a\ = a + 6, 2a s +i = c + d. 

Tratnik showed that these polynomials contain multivariable Jacobi, Meixner-Pollaczek, 
Laguerre, continuous Charlier, and Hermite polynomials as limit cases, and he used a per- 
mutation of the parameters and variables in (1.2) and (1.4) to show that the polynomials 

(1-7) 

W n (x) = W n (x; a, 6, c, d, a 2 , a 3 , . . . , a s ) 

= w ni (x 1 ;c + a 2 , s + N 2 , a , d + a 2)S + N 2 , 8 , a, b) 

s 

\ J Wn fe (^;c+a)t + i, s + iV fc+ i jS ,d + a fc+ i, s + N k+ljS ,a k + ix k -i,a k - ixk-i) 



x 



k=2 



also form a complete system of multivariable polynomials of total degree N s in the variables 
Vk = x % k = 1, . . . , s, that is orthogonal with respect to the weight function p(x) in (1.5), 
and with the normalization constant 



(1.8) 



s 

j J n k \(N kjS + N k+liS + 2a kjS+1 - l) nfc 



fc=i 



x 



r(N kyS + iV fc+M + 2a fc+M+ i)r(n fc + 2a fe ) 
r(2.V,,,.. - 2o/,,,. ,) 

x T(a + c + a 2 , s + iV s )r(a + d + a 2 , s + iV s ) 
x r(6 + c + a 2 , s + iV s )r(& + d + a 2iS + N s ). 



The Askey- Wilson polynomials denned as in [1] and [3] by 



(1.9) Pn(x\q) = p n (x;a, b, c, d\q) 

= a~ n (ab, ac, ad; q) n 4(^3 



q n ,abcdq n ,ae ,ae 
ab, ac, ad 



where x = cos6>, are a (/-analogue of the Wilson polynomials (for the definition of the 
(/-shifted factorials and the basic hypergeometric series 4 3 see [3]). These polynomials 
satisfy the orthogonality relation 



(1.10) 



Pn{x\q)p m {x\q)p(x\q)dx = X n (q)5 r . 



with max(|(/|, \a\, \b\, \c\, \d\) < 1, 



(1.11) 



p{x\q) = p(x;a,b,c,d\q) 



2i0 o -2i0. 



q) 00 ^-x 2 )- 1 / 2 



{ae id , ae~ id , be i0 , be~ ie , ce iB , ce~ i0 , de i0 , de~ ie ; q) c 



and 



(1.12) 



A n (<?) = A n (a,6,c,c%) 

27r(abcd; q) c 



(q, ab, ac, ad, be, bd, cd; q)^ 
(q, ab, ac, ad, be, bd, cd; q) n (l — abcdq~ x ) 



x 



(abedq 1 ;q) n (l — abcdq 2n x ) 



In this paper we extend Tratnik's systems of multivariable Wilson polynomials to 
systems of multivariable Askey- Wilson polynomials and consider their special cases. Some 
(/-extensions of Tratnik's [9] multivariable biorthogonal generalization of the Wilson poly- 
nomials are considered in this Proceedings [4]. (/-Extensions of Tratnik's [11] system of 
multivariable orthogonal Racah polynomials and their special cases will be considered in 
a subsequent paper. 

2. Multivariable Askey- Wilson polynomials. In terms of the Askey- Wilson poly- 
nomials a (/-analogue of the multivariable Wilson polynomials can be defined by 

(2.1) P n (x|<?) = P n (x; a, b, c, d, a 2 , a 3 , . . . , a s \q) 

"s-l 



II Pn k aA 2 , k q N «-i , bA 2 , k q N ^ , a k+1 e ie ^ , a k+1 e~ i9 ^ \q) 
fc=i 

x p ns (x s ; aA 2yS q Ns - 1 , bA 2yS q Ns ~ 1 , c, d \q) 

3 



where x k = cos/9 fc , A jyk = JJ a i; A k+1 ^ k = 1, A k = A ljk , 1 < j < k < s. Our main aim 

i=3 

this section is to show that these polynomials satisfy the orthogonality relation 

(2.2) J - J P n (x|g)P m (x|g) / o(x|g)da;i • • -dx s = X n (q) <5 n ,m 

with max(|g|, |a|, |6|, |c|, |d|, ja,2j , • • • , \ a s \) < 1, 
(2.3) 

p(x\q) = p(x; a, 6, c, d, a 2 , a 3 , . . . , a s \q) 

= {ae ie \ae- ie \be ie \be-^;q)^ 

(e 2ie « 1 e- 2ie «;q) 00 (l-x 2 k )- 1 / 2 



m 



s-l 

n 



k= 



L (a fc+ ie iefc + 1+iefe , afc+ie i6lfc + 1 _i9fe , a k +ie i0k ~ i0 k+i ? afc+ie - ie fc+i ; <?) c 



x 



(e 2 ^,e- 2 ^; 9)^(1 -x^)- 1 / 2 
(ce i6ls , ce~ l0a , de lSs , de~ l6s ; 5)00 



(2.4) 



A n (?) = A„(a, 6, c, d, a 2 , a 3 , . . . , a s |<?) 



= (2tt) 



n 



(9, A 2 k+1 q 



Nk+N^-l. 



q)n k (Al +1 q 2Nk ;q) 



(q t Alq N *+ N *-i,al +1 q»*;q) 00 
x (acA 2 , s q Ns , adA 2 , a q N ' , bcA^a^ 1 ' , bdA 2 , s q Ns ; g) ~ x , 

where a 2 = a& and a 2 +1 = cd. The two-dimensional case was considered by Koelink and 
Van der Jeugt [6], but they did not give the value of the norm. First observe that by 
(1.10) — (1.12) the integration over x\ in (2.2) can be evaluated to obtain that 



(2.5) 



J Pm(si;. 



= 6. 



a,b,a 2 e l ° 2 ,a 2 e l ° 2 \q)p mi (xr, a, b, a 2 e 1 " 2 , a 2 e ^ 2 |q) 

x p(x\;a,b, a 2 e ld2 ,a 2 e~ ld,2 \q)dx\ 

2n(q, a6aig ni_1 ; g) ni {aba 2 2 q 2ni ; 5)00 



m ,mi 



(q,abq n \alq n ^;q) oc 
x (aa 2 g n V e2 , aa^e"^ 2 , 6a 2 g n V 92 , 6a 2? ni e -^ 2 ; g)" 1 . 

After doing the integrations over xi, X2, • • • , x_j for a few j one is led to conjecture that 

(2.6) 

P^\x\q)P^(x.\q)p^(x.\q)dx 1 ■ ■■dx j 

(q, A 2 +1 q N *+ N ^-i;q) nk (A 2 +1 q 2N * ; q ) a 



= (2ny 



IP 



n k ,m k 



(?, A 2 q 



2 n N k +N k ^-L 

n.. - n n k • 



a k+ iQ nk ^)c 



x {aA 2J+1 q N i e ie ^ , aA 2ij+1 q N t e~ i9 ^ , bA 2ij+1 q N * e ie ^ , bA 2J+1 q N ^ e~ i9 ^ ; ? ) 1 , 
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where 



3 

PtfHx-lQ) = II ( x *> aA 2ik q Nk -\bA 2 , k q Nk -\a k+1 e ldk +\a k+1 e- 10 ^ \q), 
k=i 

pW(x\q) = (ae l0 \ae- l0 \be l0 \be- t01 ;q)^ 

3 i O.4.O.. 



X 



tt (e^^e-^-qUjl-xl)-^ 

(a k +ie i9k + 1+i9k , a k+ ie i9k + 1 ~ iek , a k+ ie i9k - idk + 1 , a k+ ie- i9k + 1 - i9k ; q) 



for j = 1, 2, . . . , s — 1. To prove this by induction on j, suppose that j < s — 1, multiply 
(2.6) by the Xj+i -dependent parts of the weight function and orthogonal polynomials, and 
then integrate with respect to Xj + i to get 

(2.7) 



{t, (Q,Al +1 q Nk + Nk -^-q) nk (Al + iQ 2Nk -,Q) 
lJ> fc ,™ fe - {qi Al q N k+ N k _ lia 2 + iqnk . q)oo - 



Jfe=l 
1 



x J ^ Pn J+1 (xj+i; aA 2iJ+1 Q^ , 6A 2;J+1 ^ , a J+2 e^+ 2 , a J+2 e ^+ 2 |g) 

x Pm j+1 (xj+i] aA 2 , J+1 q Nj , bA 2jJ+1 q N i , a J+2 e^+ 2 , a J+2 e-^+ 2 \q) 
x aA 2iJ+1 ^, L4 2J+1 ^, a j+2 e id >+ 2 , a J+2 e-^+ 2 \q)dx j+1 

_ (9 v+1 [ft x ( g ^ 2 fc+ i^ +jVfc - 1 - 1 ;g)n fc (Al +1 ^; 9 ) 00 1 
U " j ll d — ( q ,Alq" k +" k -^al +1 qr> k ;qU 



_k=l 



x (aA 2iJ+2 ^+ 1 e^+ 2 , aA w q N ^ e"^+ 2 , 6A 2ii+2 g Jv ^ 1 e^+ 2 , bA 2jj+2 q N: >+ 1 e"^+ 2 ; g) J , 

which is the j — > j + 1 case of (2.6), completing the induction proof. 
Now set j = s — 1 in (2.6) and use it and (2.5) to find that 

(2.8) J ---J Pn(^\q)Pm(^\q)p(^\q)dx 1 ■ ■ ■ dx 



_ r9 Y -i \ s ri, (q,Ai +1 q^ Nk -^q) n Ml +1 q 2N ^qU 

x J Pns (x s ;aA 2:S q N °-\bA 2:S q N °-\c,d\q) 
x p ms (x s ; aA 2 , s q N °-\ bA 2 , 8 q N '- 1 



■III -j • ••-"■2,s9 S S C, 

(ce i6ls , ce~ lBs , <ie i6ls , de~ t9s ; g)^ 
x (aA 2iS ^- 1 e l9 %aA 2iS ^- 1 e-^%6A 2iS ^- 1 e l9 %6A 2iS ^- 1 e-^;Q)- 1 ^ s 



x 



= A n (<?) 5n,m, 



where X n (q) is given by (2.4). This completes the proof of (2.2). 

Note that the integration region and weight function in (2.2) and (2.3) are invariant 
under the permutation of variables and parameters 

(2.9) a <-> c, b <-> d, a k+1 <-> a s - k +i, k = 1, 2, . . . , s - 1, 

0fc<->0 a _fc+i, fc = l,2,...,s. 

Hence, when these permutations are applied to (2.2) and (2.3) the transformed polynomials 
also form an orthogonal system with the same weight function. Since the polynomials 
Pn(x\q) in (2-1) are not invariant under (2.9), we obtain a second system of multivariable 
orthogonal Askey- Wilson polynomials, which is a (/-analogue of Tratnik's second system 
(1.7) of multivariable Wilson polynomials. After doing the permutation n k <-> n s - k+ i, 
k = 1, . . . , s, the transformed polynomials and the normalization constant are given by 

(2.10) P„(x|g) = P n (x.;a,b,c,d,02,a 3 ,...,a a \q) 

= p ni (xi; cA 2 , s q N2 ' s , dA 2 , s q N2 ' s , a, b\q) 

s 

II Pn fc (x fc ; cA k+lyS q N ^- , dA fc+lifl g JVfc + 1 - , ae^" 1 , ae"^" 1 |g) 



.fc=2 



(2.11) A n (g) = A n (a, 6, c, d, a 2 , a 3 , . . . , a s |<?) 



= (2t) 



n 



((/ .A^ - /7-^fe,s+-^fe + l,s" 



1. 



x (acA 2iS ^, adA 2 ,sq N ', bcA 2 , s q N % bdA 2iS q N °; g)" 1 , 

with af = a&, a 2 +1 = c<i, and max(|g|, |a|, |6|, |c|, | a-2 1 , | ci3 1 , . . . , |a s |) < 1. These polyno- 
mials are of total degree N s in the variables x\, . . . , x s and they form a complete set. 

A five-parameter system of multivariable Askey- Wilson polynomials which is associ- 
ated with a root system of type BC was introduced by Koornwinder [7] and studied with 
four of the parameters generally complex in Stokman [8] . 

3. Special Cases of (2.2). First observe that the continuous dual g-Hahn polynomial 
defined by 



(3.1) 



d n (x;a, 6, c\q) = a n (ab, ac; q) n 3 </> 2 



q ,ae , ae" 
ab, ac 



is obtained by taking d = in (1.9) and x — cos 9. Since d n (x; a, 6, c\q) is symmetric in its 
parameters by [3, (3.2.3)], we may define the multivariable dual g-Hahn polynomials by 



(3.2) 



D n (x\q) = L> n (x; a, 6, c, a 2 , a 3 , • • • , a s \q) 



s-l 



l\d nk (x k ; a k+1 e ie ^ , a k+1 e~ ie ^ , aA 2jk q N ^ \q) 



.k=l 



x d ns (x s ;b, c,aA 2 , s q Ns - 1 \q), 
6 



with Xk = cos6>fc for k = 1, 2, . . . , s. It follows from the 6 = case of (2.2) — (2.4) that the 
orthogonality relation for these polynomials is 



(3.3) 



(x.\q)D m (x.\q)p(-x.\q)dxi ■■■dx s = X n (q) 5 U)1 



with 
(3.4) 

p(x\q) = p(x; a, 6, c, a 2 , a 3 , . . . , a s |<?) 
= (ae^,ae-^; ? )- 1 



x 



X 



(e 2 *\ e- 2 ^; «,)«,(! -x^)" 1 / 2 



s-l 

n (a fc+ ie iefc + 1+iefc , a k+1 e i9 ^~ i9 ^ , a k+1 e i9 ^- i9 ^ , a^ie"^ 1 "^ ; q) c 



(be l9a ,be %9a , ce lds , ce l6,s ;q , ) o' 
(3.5) 

A n (<?) = A n (a, 6, c, a 2 , a 3 , . . . , a a \q) 



U(q nk+ \at +1 q^;q)^ 



.k=i 



(abA 2 , s q N %acA 2jS q N °;q) 00 1 , 



where a 2 +1 = be and max(|(/|, \a\, \b\, |c|, 1 0,2! , 1 0,3 1 , . . . , \a s \) < 1. 

By taking the limit a — > in (3.2) — (3.5) we can now deduce that the multivariable 
Al-Salam-Chihara polynomials defined by 

(3.6) S n (*\q) = Sn(x; 6, c, a 2 , a 3 , . . . , a s \q) 



s-l 



J\pn k (x k ; a k+1 e i9k ^ , a k+1 e~ i9 ^ \q) 



.fe=i 



x p ns (x a ;6,c|g). 
satisfy the orthogonality relation 



(3.7) 

with 
(3.8) 

p(x\q) 



J - J S n (x\q)S m (x\q)p(x\q)dx 1 ---dx s = X n (q)5 r 



(e 2iek ,e- 2i9k ;q) O0 (l-xl)- 1 / 2 



s-l 

n (a k +ie i9k + 1+i9k , afc+ie^+i _iefc , afc + ie i6,fc_i6,fe + 1 , a/ £+ ie- i9fe + 1 _iefe ; q) c 



x 



{be %9a ,be l9a , ce^ s , ce t9a ;q)oo' 



(3.9) 



A n ( g ) = (27r) s n(? nfe+1 ,4 + i9 nfc ;9)^, 



fe=i 



7 



where a 2 +1 = 6c, max(|g|, \b\, \c\, \a2\1 1 0,3 1 , . . . , \a s \) < 1, and the Al-Salam-Chihara poly- 
nomial p n (x; 6, c\q) is defined by 



(3.10) 



p n (x;b, c\q) = b n (bc;q) 



n 3<P2 



be iB M 
6c, 



see [5, (3.8.1)]. 
Setting 

(3.11) a = q^)/\ b = q^)/\ c = -g^V 4 , d = -g^+ 3 )/ 4 

in (2.1) and (2.2) gives a multivariable orthogonal extension of the continuous g-Jacobi 
polynomials Pn'^\x\q) defined in [3, (7.5.24)], while setting 

(3.12) a = q 1 / 2 , b = q« +1 / 2 , c = -q? +1 / 2 , d = -q 1 ' 2 

in (2.1) and (2.2) gives a multivariable orthogonal extension of the Pn (x; q) polynomials 
defined in (7.5.25). Also, via [3, (7.5.33)] and [3, (7.5.34) with q -> q 1 ' 2 ] the a = f3 = A-l/2 
substitution gives a multivariable orthogonal extension of the continuous g-ultraspherical 
polynomials C n (x;q x \q). By letting A — > 00 when we use (3.12), i.e. set a = —d = q 1 ! 2 
and b = c = 0, we get a multivariable orthogonal extension of the continuous g-Hermite 
polynomials defined in [3, Ex. 1.28]. 



A multivariable orthogonal extension of the continuous g-Hahn polynomials defined 



by 



(3.13) p n (cos(6 + (/>)• a, b\q) 

= (a 2 ,ab,abe 2i ^;q) n (ae^Y 



4(P3 



q- n ,a 2 b 2 q n - 1 



ae 



a , ab, abe 



2i<f>+i6 

•j 

1id> 



ae 



;q,q 



see [3, (7.5.43)], is obtained from (2.1) — (2.4) by replacing a, 6, c, d, 9k and Xk = cos 6k 
by die*^, aie - ^, a s+ ie*^, a s+ ie~^, 6k + 4> an d cos(6'fc + 0), respectively. 

It is clear that similar special cases of the second system of multivariable orthogonal 
Askey- Wilson polynomials can be obtained by appropriate specialization of the parameters 
in (2.10) and (2.11). Additional systems of multivariable orthogonal polynomials will be 
considered elsewhere. 
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